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Abstract. The paper investigates the properties of a class of resource allo- 
cation algorithms for communication networks: if a node of this network has 
x requests to transmit, then it receives a fraction of the capacity proportional 
to log(l + x), the logarithm of its current load. A fluid scaling analysis of such 
a network is presented. It is shown that the interaction of several time scales 
plays an important role in the evolution of such a system, in particular its co- 
ordinates may live on very different time and space scales. As a consequence, 
the associated stochastic processes turn out to have unusual scaling behaviors 
which give an interesting fairness property to this class of algorithms. A heavy 
traffic limit theorem for the invariant distribution is also proved. Finally, we 
present a generalization to the resource sharing algorithm for which the log 
function is replaced by an increasing function. 
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1. Introduction 

The resource allocation problem considered in this paper involves J nodes which 
have access to a common shared resource, for example a communication channel 
or a processing unit. The latter is assumed to have a fixed capacity, say 1. The 
resource is shared among nodes in the following way: for 1 < j < J, if node j has 
rij requests pending, it receives the instantaneous fraction of capacity 

(i) IM 

from the resource. The algorithm is thus defined by the function x t— > f(x). There 
are several situations where the capacity is allocated in this way. 
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1.1. Saturated Node of the Internet. In this context, the nodes correspond 
to TCP flows with different sources and destinations. The resource here is the 
processing time of a fixed router on the path of these flows. The packets of a flow 
are queued in the buffer of the router until they are routed to the next stage of 
their path. Congestion is simply the situation when the buffer is full and therefore 
incoming packets are lost. Because of the TCP protocol, a given flow will increase or 
decrease the rate at which it sends the packets, depending on the level of congestion 
of the routers on its path. There are several ways to represent this phenomenon. 
It should be kept in mind that the following descriptions are mathematical models 
of the way TCP is thought to allocate the bandwidth, not of the TCP algorithm 
itself. See Massoulie and Roberts [15] . 
a) Processor- Sharing disciplines. 

A popular, simplified, stochastic model of this situation consists in consid- 
ering that the router allocates its processing power to each flow according 
to a slight generalization of the allocation policy given by Relation (JlJ with 
a function / depending on the node j and of the form Wjn, where 1/iVj can 
be the round trip time between the source and the destination. This alloca- 
tion algorithm corresponds to the Discriminatory Processor-Sharing Policy. 
Node j has an instantaneous fraction of capacity given by 

w jnj 



(2) 



WiUi + w 2 n 2 + ■ ■ ■ + wjnj 

See Altman et al. [5] and references therein for a survey. When all the 
Wj's are 1, we obtain the classical processor-sharing policy: node j receives 
the fraction of capacity nj/{n\ + • • ■ + n.j), and the bandwidth is equally 
divided among the current requests. In the last ten years, different classes 
of stochastic models of processor-sharing policies have been extensively used 
to describe the congestion in IP networks. 

b) Alpha-fair disciplines. 

These policies have also been introduced to describe the allocation of band- 
width in IP networks (see Mo and Walrand [IS]), in terms of an optimization 
problem (cf. Kelly et al. [T3] )• In our context, a related policy would cor- 
respond to the case f(n) = n a , n > 0, so that a non-empty node j has an 
instantaneous fraction of capacity given by 



n 7 



(3) 

n i 1" "2 i" • ■ ■ -r lo j 

The case a = 1 is the processor-sharing discipline presented above. 

In the wireless section below, the situation is quite different since the bandwidth 
allocation algorithm is defined explicitly by relations similar to ([I}. 

1.2. Wireless Networks. This is again a simplified, but meaningful stochastic 
model of bandwidth allocation, this time in wireless networks. The resource here 
is a radio channel in a region where there are J stations/mobiles. At a given time, 
because of interferences, only one station can transmit successfully in this region. 
A station with nj messages waiting for transmission can detect if there is a commu- 
nication going on or not. If not, a classical backoff mechanism is used: the station 
starts transmitting after an exponentially distributed amount of time with param- 
eter f(rij). If another station starts a transmission before that time, the attempt 
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of transmission is canceled. Consequently, if initially there is no transmission, node 

j will access the channel successfully with probability f(nj)/(f(m) H + f(nj))- 

See Abramson pQ and Metcalf and Boggs [19 for historical references, and Tassiulas 
and Ephremides |27j . If a small quantum 5 is transmitted at each access, it is not 
difficult to see that, provided that backoff times are small (which is the case if one 
of the components is large) , as 8 goes to the effective capacity allocated to station 
j is indeed given by 

(4) IM 



f{ni) + ■■■ + f(nj) 



This is the analogue of the approximation of the round robin policy by the processor- 
sharing discipline. 

Fair Access to Resource: The Choice of the function x i-> f(x). 

The function / should clearly be increasing, so that the fraction of the capacity 
allocated may grow with the number of requests. This is the case if f{x) = x a 
which corresponds to the Alpha-fair disciplines already mentioned. However, these 
policies may have a serious drawback. Indeed, if a station j has a large number of 
requests pending while the other stations are lightly loaded, the latter will receive 
a negligible fraction of the bandwidth. The station j will therefore capture the 
channel for its own benefit, until the instant when some of the other stations reach a 
comparable level of congestion. This is a highly undesirable property for a network 
where fairness issues (for nodes, not requests) are of primary importance. See 
Bonald and Massoulie [5]. 

A possible way of solving this problem is to consider increasing functions / which 
grow slowly to infinity like, for example, the concave function x i-> log(l + x), or 
x i Y log log(e + x). In this way, one can expect to reduce significantly the impact of 
saturated nodes even if they still receive a sizable fraction of the available capacity. 
Related algorithms have been considered in the context of wireless networks, see 
Shah and Wischik [55] and references therein. Bouman et al. [B] and Ghaderi 
et al. [T2] investigate the impact of the "agressivity" of the function / on the 
stability and on the delays for several wireless network architectures with a related 
bandwidth allocation scheme. 

In this paper, we mainly investigate the case f(x) = log(l + x). The general case 
is sketched in Section El The instantaneous fraction of capacity of the j th node is 
therefore given by 

(5) log(l + "j) 

log(l + ni)+log(l + n 2 )+---+log(l + ?v)' 

The log function moderates the rate at which a saturated station tries to access the 
resource, which is a desirable property in an heterogeneous network where traffics 
may have very different characteristics. In the context of wireless networks, a related 
algorithm was used to show that an optimal stability region is possible in a quite 
general network. The growth properties of the log function play an important role 
in the proof of the result. Basically, the log of the states of the saturated stations 
being quite stable on some large time intervals, the schedule (the set of stations that 
can transmit at some time) quickly reaches some equilibrium and stays around it. 
A Lyapounov function argument can then be used to prove ergodicity (see Shah and 
Shin [25]). Up to now, apart from these stability results, little is known about the 
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quantitative and qualitative properties of these algorithms. As we shall see below, 
the mathematical analysis of this class of algorithms presents some challenging and 
unusual problems (see also Wischik 28 ). We first achieve a fluid limit scaling 
analysis, which gives a very precise description of the qualitative behavior of these 
algorithms. Additionally, we derive a heavy traffic limit theorem result for the 
invariant distribution of the associated Markov process. Before presenting our 
main results, we briefly recall the main definitions of the fluid limit scaling. The 
interested reader will find an extended presentation in Bramson [8 or in Chapter 8 
of Robert [13]. 

Throughout the paper it will be assumed that, for every 1 < j < J, the requests 
arriving at the jth node form a Poisson process with rate Xj > 0. Each request at 
node j leaves the network when it has received an exponentially distributed amount 
of time with parameter \ij from the common resource. The average load of the jth 
node is denoted by pj = Xj//ij. 

Fluid Limits. The fluid limit scaling of a stochastic process (Z(t)) in K J consists 
in speeding up time and space in proportion to the norm of its initial state: 

Z*(f) = ^p, with TV = ||Z(0)||. 

A possible limit in distribution of the sequence of processes (Z]y(t)) is called as a 
fluid limit of the process (Z(t)). Hence, in some sense fluid limits give a first order 
description of (Z(t)). This is a convenient tool to investigate multi-dimensional 
processes for which general results are scarce. In the context of Markov processes, 
there is an additional interest since the ergodicity of the process can be connected 
to the fact that fluid limits (whose initial states lie on the unit sphere) return to 
the origin. See Rybko and Stolyar [53] and Dai [5]. Note however that the fluid 
limit scaling is well suited for processes that behave locally like random walks. For 
other processes, different scalings may have to be considered. 

In the case of the generalized processor-sharing policy defined by Relation ([2]) , 
for every 1 < j < J and t > 0, let Xj(t) denote the number of jobs waiting at the 
jth node. The evolution of this process can be represented as 

Xi (t) = A, (0) + Ma (t) + \A- u 7 / — w j*j( s ) ds 

jW 3K ' 3 ' 3 3 Jo w 1 X 1 (s)+W2X 2 {s) + ---+wjXj(s) 

where (Mj(t)) is a martingale. The scaled process with N = |A(0)| is thus given 
by 



?(t) = zf(o ) + *SM 



+ \jt-lij ZZN =A ds - 

Jo w 1 X 1 (s) + w 2 X 2 (s) H h WjXj (s) 

A standard tightness criterion and the fact that the martingale 

{(Mj(Nt)/N),l<j<J) 

converges in distribution to imply that any fluid limit ((xj(t)), 1 < j < J) should 
satisfy the ordinary differential equations 

(7) *p.ft) = A, - fjtj rT , 1< j < J, 
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on the interval [0,io], provided that the vector (xj(t),l < j < J) does not hit 
before to. See Ben Tahar and Jean-Marie [3], Ramanan and Reiman [22], and 
references therein. 

Similarly, for alpha- fair disciplines the corresponding fluid model (xj(t)) is the 
solution to the ODE 

Observe that for these two choices of function / the scaled process satisfies an 
autonomous ODE, like ©, with a stochastic noise component that vanishes as A 
gets large. Secondly, a remarkable feature of these convergence results is that all 
coordinates of the scaled process are of order A. That is, as long as x[t) is not the 
vector 0, one has Xj(t) > for all 1 < j < J. 



Problem of Fluid Limits for Algorithms with Logarithmic Weights. Let 

(Lj(t), 1 < j < J) denote the Markov process associated to the policy with loga- 
rithmic weights, i.e. associated to Relation ([5]). Due to the log function, the scaled 
process (Lj(t)) does not have an autonomous representation analogous to ©. In 
fact, it is easy to see that the corresponding stochastic equations involve both 
Lj(Nt) and log(l + Lj(Nt)), two quantities which evolve on very different scales. 
For this reason, there is no way of guessing a system of plausible "fluid equations" 
corresponding to System ([7]) [resp. to ([8])] for discriminatory processor-sharing pol- 
icy [resp. Alpha- fair policies]. See Wischik [55] and Ghaderi et al. [15]. Note that 
the question of stability of the system is not an issue here. Indeed, because of the 
work conserving property of these policies, a necessary and sufficient condition for 
the ergodicity of {Lj(t), 1 < j < J) is simply given by 

A 7 - 

Pi H h pj < 1, with pj = 1 < j < J. 

To the best of our knowledge, there is no explicit expression known for the invariant 
distribution. The fluid scaling gives a first order description of the behavior of 
this policy. As we shall see, an interesting convergence result for the invariant 
distribution just below saturation can also be derived from these results. 



Outline of the Paper. Section [2] presents the main results of the paper. Section[3] 
introduces the notation and the stochastic differential equations associated to the 
Markov process (Lj(t),l < j < J). Sections 0] [S] and H] are respectively devoted 
to the scaling properties of the time scales t h-> N f , t i-> A" 1 log At and t H> Nt. 
There we provide a precise description of the evolution of the network, together with 
some estimates of hitting times. The key results on the fluid limits are presented 
in Section [6] In Section [7] we prove a heavy traffic limit theorem for the invariant 
distribution. The case of a two node network with a general / is discussed in 
Section [5J The corresponding time scales are identified in this case. Section [9] gives 
a brief sketch of the case of a network of J nodes. 



Acknowledgments. The first author would like to thank Damon Wischik, whose 
presentation at the ICMS workshop in Edinburgh in 2010 is one of the motivations 
at the origin of this work. 
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2. Presentation of the Results 



The General Picture. 

The main result of the paper is that the states of the nodes may live on different 
space and time scales depending on the set of parameters (Xj,fij). Indeed, there 
exists a set of real numbers (ctj, 1 < j < J) with values in (0,1] such that, if 
ctj < 1, then on the fluid scale Lj is of the order of ||L(0)|| Qj and the process 
t i — ^ Lj(||L(0)|| Qj ' log ||L(0)||t) has a local equilibrium which can be expressed in 
terms of an Ornstein-Uhlenbeck process. In the next paragraph we give a more 
precise description of this phenomenon in the case J — 2. An illustration of the 
different scales is given in Figure [TJ where the y-axis is on a log -/log TV scale. 

The case of Two Nodes. 

Even in the case where there are only two classes of customers, J = 2, giving an 
asymptotic picture of this queueing system is already a challenging problem. To 
concentrate on the interesting case, let us assume that the parameters satisfy the 
following conditions: p\ < 1/2 and p2 > 1/2 and the initial state of the process is 
(L?(0),L?{0)) = (0,N). 

a) The time scale t — > N . 

A convergence result, Proposition [21 shows the equivalence in distribution 



(Lf(7V t ),,0<i<^) ~ 



Xi-Vi^-)N\0<t<al) fora^-^- 
t+l J J 1 - Pi 



The process stays at N on this time scale. The condition p\ < 1/2 
implies that a* < 1. Note that the pre-factor of iV* vanishes at t = a\. 
For this reason, this convergence result does not prove that values of the 
order of N" 1 can be reached. In fact, an extra (log TV) factor is required 
and Proposition [3] shows that if 8 < 1, the average hitting time of the value 
lSN a i\ by (if (f)) is bounded by KiN a "i log TV for some K x > 0. On the 
other hand, reaching the value [TV" 1 J is slightly longer: the average hitting 
time of [N a i\ is bounded by CiV Q i(logiV) 2 loglog(A^) for some C > 0, see 
Relation ([2SJ in Section [5j 

b) The time scale t -> N a *i (logiV) t. 

We now assume that Lf(0) = |A^J and L 2 (0) = N. Theorem [1] then 
proves the convergence of the processes 

/ L?(N°H\ogN)t)-N^ 
V y/N a 'i \ogN 

to an Ornstein-Uhlenbeck process. In other words, on this time scale is 
stabilized around the value N Ul . Again, the process stays at N on this 
time scale. 

c) The fluid time scale t —¥ Nt. 

This time scale is natural, since the network empties on it. Theorem [3] shows 
the convergence, 
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for the convergence in distribution of processes, with 

7(<) = (1 + (A 2 - M2 (l- Pl ))i) + . 

Consequently, as long as the fluid limit of (L 2 {tj) is not 0, the process 
lives on the space scale N ai . 




Figure 1. A first order picture of the network with p\ + p 2 < 1, 
pi < 1/2 and (L^O), L 2 (Q)) = 

Properties of algorithms with logarithmic weights. 

The bandwidth allocation algorithm with logarithmic weights exhibits an interest- 
ing property. For example, when J = 2 and the system is overloaded (pi + p 2 > 1), 
if pi < 1/2 the size of the queue of class 1 requests grows at a rate proportional to 
with a\ < 1. This implies that queue 1 is stable at the fluid level, i.e. that 
Li(t)/t goes to in distribution as t becomes large. Hence, without any priority 
mechanism among nodes, if a node has a light load, p\ < 1/2, then most of its 
messages will be transmitted with success even in the case where the system is 
globally saturated. Recall that in the transient case of the processor-sharing policy 
or even with the Alpha-fair disciplines, this is not true at all: the states of the 
nodes diverge to infinity at the same speed, linearly in time. This is an interesting 
feature from the point of view of fairness issues. It can be shown that an analogous 
property is valid for the network with J nodes, see Section [9] 

This study also suggests that, for general networks using such algorithms, several 
space and time scales should play an important role in the dynamics of the system. 
For the two node network described above, when the initial state is (0, N) we prove 
that the fluid limit is given by 

((0,1 + (A 2 - M2 (l- Pl ))t) + ) . 

This shows that node 2 receives the capacity 1 — p\ , which is another way of saying 
that node 1 is stable at the fluid level. The simplicity of this expression somewhat 
hides the complexity of the situation, since the quantity a* does not show up. 
Yet, as we have seen a\ has a crucial impact in the transient well as in 

the equilibrium situation: if the pair of variables (Li >p ,L2, p ) has the equilibrium 
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distribution of the Markov process (Li(t),L2(t)) and if pi < 1/2 is fixed, the heavy 
traffic limit-Theorem |4] shows the convergence in distribution 



where X is an exponential random variable. Hence at equilibrium and in the heavy 
traffic regime, the relation L\ ~ Lj 1 a l so holds. 

It is unlikely that a standard fluid analysis, i.e. deriving directly some equations 
similar to Relation for example, can be done to investigate the qualitative 
behavior of more complex networks. This is where the consideration of the various 
time scales is useful. It gives a tool to explain, via a dynamic picture, the multiple 
orders of magnitude of the state variables at equilibrium. 

An Interaction of Time Scales. 

There is an unconventional property for the fluid scaling of a queueing system. For 
most of the queueing networks investigated up to now, the classic general scheme 
for the fluid scaling of the associated Markov process (X(t)) is as follows: there is 
a subset of the coordinates whose values are of the order of N = ||X(0)|| and the 
other coordinates form an ergodic Markov process whose invariant distribution de- 
termines the evolution of the large coordinates on the fluid scale. See Malyshev [17] 
and Bramson [8] for some examples. 

Here the situation is different. If the initial state is (0, AT) and if p\ < 1/2, 
then the large coordinate Li is of the order of iV but L\ is an order of magnitude 
smaller, N ai with < a* < 1. There is indeed an underlying ergodic Markov 
process but at the second order, namely an Ornstein-Uhlenbeck process scaled by 
a factor y/N a i log TV. See Relation ©. 

The associated stochastic model exhibits a stochastic averaging principle at the 
origin of the second expansion in Relation (flO|) . The key technical result of the 
paper, Theorem[2j states that when pi < 1/2, on the fluid time scale t H> Nt, L 1 ^ is 
uniformly of the order of (L^)" 1 on any finite time interval with high probability. 
Recall that 

a) If if (0) = N and Lf (0) = N a * , then on the time scale t i-» N a i (log N)t 
we have - L$ (0) and Lf is of the order of (i^(0)) Q i. Additionally, 

can be represented by an Ornstein-Uhlenbeck process around A~ ai , see 
point b) in the previous paragraph. 

b) On the fluid time scale, (Nt) is of the order of 7(^X2 (0) with j(t) defined 
above by Relation (|T0| . 

The problem lies in proving that on the fluid time scale adapts sufficiently 
quickly to preserve the relation ~ ) ai . A central limit result, Proposition[5j 
suggests that this is not the case on the timescale of the Ornstein-Uhlenbeck process, 
at least for a second order description. On the other hand, on the fluid time 
scale Theorem [2] shows that this separation of time scales holds. Its proof uses 
several estimates related to average hitting times of reflected random walks and 
some coupling arguments. One of the problems encountered is that the potential 
natural stochastic fluctuations of the fluid time scale, of the order of y/N, can be 
large compared to A^" 1 (if a* < 1/2 for instance). In particular, standard stochastic 
calculus cannot be used as such to prove the result. It turns out that the potentially 
large fluctuations are reduced by the strong ergodicity properties of the underlying 
Ornstein-Uhlenbeck process. Thus, it does not seem that the classical techniques to 
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prove stochastic averaging results can be used here. See Khasminskii [14] , Freidlin 
and Wentzell [H] and Papanicolau et al. [21] for a general presentation of methods 
to prove stochastic averaging principles. 

3. The Stochastic Model 

In this section, we introduce the main stochastic processes and some notation. 
If h is a nonnegative Borelian function on R + , we let Nh denote a Poisson process 
with rate x H > h(x) on R + . This process can be defined as follows. If V is a 
homogeneous Poisson point process on with rate 1 and / is some nonnegative 
Borelian function on R + , then Afh(f) is defined by 

f(u)Af h(u) (du) = f f(u)V([0, h(u)] x du). 

Jr 2 + 

For £ > 0, Af% denotes the Poisson process with rate £ on R + , i.e. corresponding to 
the constant function equal to £. In addition, for any < a < b, A/^([a, 6]) stands 
for the number of points of A/j in the interval [a, b\. Throughout the paper, the 
various Poisson processes used will be assumed to be independent. 

We consider two classes of customers. The arrival process of class j customers 
is a Poisson process with rate Aj, the distribution of the duration of the required 
service is exponential with rate fj,j, and pj denotes the ratio Xj/fij. Each class of 
customers has a dedicated queue and there is a single server working at unit speed. 
If the state of the system is (x\, x^) € N 2 , where Xj is the number of jobs in queue 
j, then customers of class j receive the fraction of service 

M-n wr ^ dcf. log(l + z. t ) . 

11 Wi{x 1 ,x 2 ) = : — 7— r— j — 7T- c> « = 1,2, 

log(l + xx) +log(l + x 2 ) 

from the server, with the convention that 0/0 is 0. The process of the number 
of jobs in queue j G {1,2} is denoted by (Lj(t)). Since we are only interested in 
the total number of customers of each class, there is no need to specify the service 
discipline for each queue. It can be Processor-Sharing or FIFO (First In First Out), 
for example. 

Stochastic Differential Equation. The stochastic process (Li(t), L2(t)) can be 
expressed as the solution to the following stochastic differential equation (SDE): 

(12) dL i (t)=7V A< (dt)-^ jm(il(t _ ) ,i 2 ( i _ ) )(dt), i = l,2, 

where Li(t—) denotes the left limit of Li at t and Wi is the function defined by 
Relation (fTT) . 

A Saturated System. For N 6 N, A, /i > 0, it will be convenient to introduce 
a one-dimensional process (Xpf(t)) describing the evolution of the number of cus- 
tomers in a given queue when the number of jobs in the other queue is "large", i.e., 
of the order of N. The process (X/v(i)) is thus defined as the solution to the SDE 

(13) dX N (t) = Af x (dt) -^ Wl{XtN _ 1} (dt). 

From a Markov process point of view, (Xjv(t)) is simply a birth and death process 
on N whose Q-matrix (Q(x, y)) is defined by 

q(x, x + 1) = A 

/ -, \ log(l + x) 
q(x,x — l) =u-, — T- ^ — ; 777, x > 0. 
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As we shall see, when (£i(0), L 2 (0)) = (Q,N- 1), X N (0) = 0, A = Ai, and p = px, 
the two processes and (XN{t)) are close enough (for our purposes). Note 

that this is not completely clear since the process (L 2 (t)) may drift away from N 
and therefore change the service rate received by each class. It turns out that, 
because of the slow increase of the log function, this property will hold at least at 
the beginning of the sample paths. 

By integrating the SDE (fT3|) one obtains that, for any t > 0, 

(14) X N (t) =X N (0)+M x ([0,t])- f ^ Wl (x N (u-),N-i)(du) 

Jo 

= X »W + ^ - p f . ,!f * + ^ + ] N du + M »®> 
J log(l + X N (u)) + logTV 

where (Mjv(i)) is the martingale 

M N (t)=M x ([0,t])-Xt + f [^ l( x JV («-),Jv-i)(du) - fiWx(X N (u),N-l)du], 

Jo 

whose increasing process is given by 

(15) (M N ) (t) = Xt + p — - ... — — du. 



I log(l + X N (u))+logN 
4. The Initial Phase 

This section is devoted to the very beginning of the evolution of the first com- 
ponent (Lx (t)), when it starts from while (0) = TV. To start with, we have 
the following asymptotic result on the initial growth rate of the process (X^(t)) 
defined by Equation (|13[) . Here and later, we write a A 6 for the quantity min(a, b). 

Proposition 1. If Xn(0) = and 

* def. p . A 

a = , where p = — , 

1-p p 

then the convergence in distribution of stochastic processes 

li m ( Xn ^ J <t <a* Al) = (x- pt-^—, < t < a* A 1 
W^ + oo \ iV* / v t + 1 

ZioZds for the uniform topology on compact sets of (0, a* A 1). 



Proof. The evolution equation (|14j) and a change of variables give us that for every 
t > 0, 

XNjN*) = X N (1)-X-M N (1) M N {N*) 
N* N* N* 

A J \og(l+X N (N"))+logN K S 7 
Letting Z N (t) d = (1 + ^(iV 4 ))/^*, we thus have 

(i6) Mt) = ^)iy-^) + w 

+ X — p f ^(Z N (t V ))Ht- V )logN 

' io log(Zjv(t -«)) + (<-« + !) log JV V 8 j 
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Let us first show that the martingale term does not play a role for this scaling. 
Indeed, for < a < b < 1, Doob's Inequality yields for every e > 0, 

SU pM >£ )< P ( sup \M^M >£ ) 

< ^Ta E (\MN(N b )\ 2 ) = ^MM N ) (N b )) < ^±^N b ~ 2a . 

The last term can be made arbitrarily small when N is large by choosing b < 2a. 
Since the interval (0,a*) can be covered by a finite number of such intervals, the 
martingale term is indeed negligible with probability tending to 1 as N goes to 
infinity. The relation Xjy(l) < Af\([0, 1]) implies that the first term in the r.h.s. of 
(fTBj) vanishes too when divided by AT*. 

Next, the inequality Xjv(s) < Af\([0, s]) gives us that 

\og(Z N (t-v)) + (t-v + l)\ogN [ ' 0gly)ly 
, v fA dcf. f* log((l+AA A ([0, N^/N^+jt-v) log N 
~ YN{t) = Jo log((l+AA A ([0,iV-]))/iV-') + (^+l)logiV (1OgiV)iV ^ 

By the law of large numbers for Poisson processes, for any < a < b the process 
(Af\([0,N s ])/N s , a < s < b) converges in distribution to (A, a < s < b). As a 
consequence, 

lim (Y N (t), a<t<b) = ( : -, a<s<b 

N^+oc \t + 1 

Gathering these estimates, we obtain that for any < a < 6 < a* A 1 and any 
e > 0, 

lim F\\-li-^ e< inf Z N (s) < sup Z N (s) < A + z | = 1. 

iV-f+oo Y 0+1 a<s<b a<s<b J 

Since A — + 1) > for 6 < a*, this shows that the term log(Zjv(i — v)) in 

the integral appearing in the r.h.s. of Relation (|16p can be neglected in the limit 
on any interval ue [0, t — a]. But on [t — a,t], the term multiplying (log N)N~ V is 
bounded by 1 and so we obtain the desired convergence in distribution. □ 

The basic ingredient of the following proposition is that, under the initial con- 
ditions (if (0),if (0)) = (0,A0 and X N (0) = 0, then the processes (if (i)) and 
(X (t)) are close on the time scale t i— > AT', < i < a* and therefore that the 
convergence result of Proposition [1] is also valid for the process (if (t)). 

Proposition 2. //(if (t),if (t)) is i/ie solution of the SDE (JT2J) wit/i initial con- 
dition (0, AT) and a 1 =pi/(l — pi) urai/i pi = Xi/jix, then the convergence 

fL N (N v ) \ { t 
lim 1 y t ' , 0< i < a? A 1 = Ai - pi , < t < at A 1 

N^ + oo \ AT* 1 / V t + 1 1 

/io/ds /or t/ie uniform topology on compact sets of (0, a* A 1). 

Proof. Since if (0) = A" and the number of jobs of class 2 decreases at rate at most 
\i\ and increases at rate Ai, for any < b < 1 there exist two constants < rj < 7 
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such that if 

An \ VN < inf L% (a) < sup LF(a) < ~/N I , 

^ 0<s<N» 0<s<N" J 

then P(An) tends to 1 as N tends to infinity. On the set An, the jump rate for 
departures of lies between fi{W^{-,N) and fnWi(-,N), where 

Wf{x,N)- 1 ° S{1 + X) 



io g (i + x) + iog(<yjv) 

Now if (X N (t)) denotes the solution to Equation (fl"3|) with W\ replaced by Wf, a 
straightforward coupling shows that on the set An the relation 

X N (s) < Lf (a) < X N (s), 0<a<N b 

holds almost surely. A glance at the proof of the convergence result of Proposition!]] 
shows that this result also holds for both processes (X N (s)) and (X N (a)), and so 
the proposition is proved. □ 

The above proposition shows that if a* < 1 (i.e. p\ < 1/2), then on the time 
scale 1 1-)- N* , < t < a* , we have 

Lf(iV*)~ (Ai-Mi^) NK 

In particular, the process reaches the quantity N ai ~ £ for any < e < a*. Note 
that when t /*a* the quantity multiplying N l vanishes, so that this convergence 
result does not show that the value iV a i is indeed reached. In Sections [5] and |51 
we shall prove that the process lives in fact in a "small" neighborhood of N Ul . 
This local equilibrium around TV" 1 is the key phenomenon to grasp in order to 
understand this bandwidth sharing policy. For now, we conclude this section by 
proving that for any < 8 < 1, the value 5iV a i A1 is attained. This is done by 
providing an estimation of the corresponding hitting time. 

Proposition 3. With the same notation and assumptions as in Proposition^ and 
if H a denotes the hitting time of a > by L^, 

H a = mf{t > : Lf (t) > a], 

then 

a) if a* < 1, for any < 5 < 1 there exists a constant C > such that for any 
N > 1, 

b) if a\> 1, for 6 > sufficiently small we have 

limsup ^ E {H 5N ) < t — -77: ■ 

Proof. As we did for Proposition^ let us prove these two inequalities for the process 
(-Xjv(t)). We use the simplified notation of Proposition [1] 
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(17) 



Let us first assume that a* < 1. For x > 0, the elementary relation 
log(l + x) a* 



log(l + x) + log A a* + l 

log((l+x)/A Q *) 



{a* + l)(logA)(l + a* + log((l + :r)/A Q *)/(logA)) 1 

together with the identity 

x = Mr- — r 

1 + a* 

and Equation (ITU) written at time iV Q (log A) i give the representation 
(18) Z N (t) d = ^ Na '0°BN)t) _ M N (N a ' (log A)i) 



log(iV- a * +Zat( U )) 



• du. 



(1 + a*) 7 a* + 1 + log(7V- Q * + Zjv(u))/(log 
Let 

rjv = inf {s > : A w (A a * (log A>) > <5A Q *}. 
From Doob's optional stopping time theorem and the fact that Zn(tn At) < 
A~ Q [(5 TV" ] , we obtain the inequality 

r^ a 'i , p iog(g+2Jv-°*) 

A«* (l + a*)a* + l + log(<5 + 2A-«*)/(logA) 1 ^ " 

Since H SNa * < N a (log N)tn, item a) is proved. 

Assume now that a* > 1, i.e. that p > 1/2. Equation (|14l) written at time Nt 
gives the relation 

" Nt log(l + A0v( S )) 



X N {Nt) = M N (Nt) + XNt - fi 
from which we can write that 



o log(l+Ajv( s ))+logA ' 

io g (i+r<wi) 



[6N] > E(X N (H SN A (Nt))) > E(H 5N A (Nt)) yX ^ + m)) + ^ 
Letting t go to infinity, the monotone convergence theorem gives us that 
(19) ]i m8U1? lE(H SN )<—^—. 

If we now choose S sufficiently small so that, with high probability, the component 
L2 (t) is still of the order of N at time 6N/(X — /z/2), a straightforward coupling 
between and some X% (recall the notation X 1 ^ from the proof of Proposition [2]) 
ensures that (fTi?)) holds as well for the process (L^(t)). This completes the proof 
of Proposition [3] □ 

5. A Local Equilibrium 

This section is essentially devoted to the behavior of our two-dimensional process 
on the time scale t n- A Q i(logA)i, when the initial state is (if (0), £f(0)) = 
(<5A Ql , A) and p\ < 1/2. The following result, Proposition [4j shows that on this 
time scale the sample paths of (L± (t)) have values of the order of xN a ^ , where < 
x < 1. When the initial value of (L^(t)) is A" 1 , we shall prove in Theorem [1] that 
the process is stabilized around A Ql . At first sight, the interest of Proposition!!] and 



14 



PHILIPPE ROBERT AND AMANDINE VEBER 



of the associated central limit theorem (Proposition [5] below) may seem marginal. 
This is not true at all since, as we shall see in Section HI the process is also of the 
order of ^(t)N a ^ on the fluid time scale 1 1— > Nt. Furthermore, the main difficulty 
in the key technical result of Theorem [2] is precisely connected to the interaction of 
these two time scales t i-» iV" 1 (log N)t and t h > Nt. 
Recall the notation a* = pi/(l — pi). 

Proposition 4. // p x < 1/2 and (L^(t), L% (t)) is the solution to the SDE (0 
with initial conditions L 2 (0) — N and L 1 ^ (0) ~ SN a i for some 8 £ (0, 1], then the 
sequence of stochastic processes 

Lf (N< (log N)t) 
converges in distribution to (h(t)) defined by 



(20) 



h=l if <5 = 1, 



h(t) 



du = ~ n ^ *\? l f 5 ^ L 



log(u) (1 + al) 



Proof. As we did for Proposition [2j the convergence is proved for the process 
(Xj^(t)). The result for (L^(t)) follows from a similar coupling argument. 

Let us first show that if X N (0) = [5N a * \ , then X N /N a * remains within [5/3, 3] 
on a given time interval [0, N a (log N)T] with probability tending to 1. Indeed, 
writing H a for the hitting time of a by X/v(-/V Q (logiV)-), the strong Markov prop- 
erty of Xn gives us that 

(21) P(H ma * <T) =E(l {H2Na . <T} F l2Na . ] (H 3Na * <T-H 2Na *)) 
<E(l {ff2NQ , <T} P r27VQ . 1 (i7 3 , VQ , <T)). 

Now, due to the monotonicity properties of the service rate of Xn , if Xn (0) = 2 N a 
then we can couple Xn with the process (2N a + i?/v(t)) defined by 

— Rn Rn + 1 at rate A, 

— Rn — > Rn — 1 at rate 

log(27V Q *) a* C 



log(27V Q * ) + log N ^ a* + 1 log N ' 

for some C > 0, 
— Rn(0) — and Rn is reflected at 0, 

in such a way that XN(t) < 2N a + i?Ar(i) for every t > 0. Hence, there remains 
to prove that Rn does not reach N a in less than N a (log N)T units of time. But 
by Kingman's inequality (see Kingman [15 or Relation (3.3) of Theorem 3.5 in 
Robert [23]) and the fact that A = pa* / (1 + a*), if On stands for the time span of 
the first excursion of Rn away from we have 

C'N°' 



sup R N (s) > N a < exp 
se[o,e N ] J \ logN 

for some explicit C" > 0. Since the ith pair of consecutive excursions is separated 
by an Exp(X) amount of time Ei and since these exponential times are independent 
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of each other, we have 

P( more than 2AiV Q * (log N)T excursions before time N a " (log N)T ) 

(2AiV a * (log jV)T \ 
J2 Ei< N a " (log N)T\ < e - c " Na * ^ N ^ T 

for some C" > 0. Coming back to (|2~Tj) . we obtain that 

P(# 3JV «* <T) < e - c " WQ *( logAr ) T + 2A7V Q *(logiV)rexp|-^^|, 

which tends to as N tends to infinity. Finally, since the infinitesimal drift of Xm is 
positive when Xn < N a , the same method can be used to show that Xn remains 
above (6/ 3) N a on the interval [0,N a (logiV)T], with overwhelming probability. 

Let now (w/(£)) denote the modulus of continuity of the function (/(f)) on [0, T}. 
That is, 

w f (0= sup |/(s)-/(t)|. 

s,t<T 
\a-t\<£ 

Let us also write again 

dcf. X N (N a * (log N)t) 

N{ ) jy a . 

Using the bounds on Xn we have just obtained, we can deduce the existence of a 
constant A independent of N such that with probability tending to 1 we have for 
any s <t E [0,T], 

log(iV- Q * + Z N (u)) 



du < A(t - s). 



a* + 1 + log(iV-«* + Z N (u))/(logN) 

Together with Relation (TT5)) and the fact that its martingale term vanishes as N 
gets large, we can conclude that for any e > and rj > 0, there exists £ > such 
that 

P(u>z*(0 >V)<£ 

holds for all N. By the tightness criterion of the modulus of continuity (see Theorem 
8.3 in Billingsley [4]), the sequence of processes {Z^{t)) is thus tight and any 
limiting point h satisfies the relation 

h(t)=S- '\ [ log(%))du. 
(1 + a*) z J 

Proposition 0] is proved. □ 
The above proposition can be seen as a kind of law of large numbers, with 
Lf(7V Q i log(iV)t) - h(t)N a "i , as N -> oo. 
It is thus natural to expect for the corresponding central limit theorem that 

/ L?(N a i log(JV)t) - h{t)N a i \ 
V ^jN a l \og(N) J 

should converge in distribution to some Gaussian process (R(t)). However, the 
following proposition shows that such a convergence cannot hold: the centering 
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term h(t) has to be replaced by a deterministic term /ijv(i) which deps on N. As 
we shall see, /ijv is such that the following expansion holds 

(a* + l) d logA' J \logN 

In particular, (/ijv(i)) converges "slowly" to (h(t)) at the rate 1/ log AT instead of a 
rate l/N ai ^ 2+e for which a "classic" centering procedure would be enough. In the 
end, this second limit theorem gives the representation 



Li(iV Q i log(JV)i) ~ h N (t)N a * + y/N^UogN R(t). 

Proposition 5 (Central Limit Theorem). If pi < 1/2 and (L^(t), (i)) * s ^ e 
solution to the SDE (|12j) wi£/i £/ie initial conditions (0) = -/V and (0) being 
such that 

if (0) - 6N a " 
am — — = y, 

TV^+oo ^iV"^ log TV 

/or some < 5 < 1 and t/elR. Then, we have 

hm f^(^do g ^)-^(^^ m))i 

/or f/ie convergence in distribution, where (/ijv(i)) is i/ie solution to the ordinary 
differential equation 

; /.x lograjv(i) 



(1+aJ) K + l + log(M*))/logA0' 
urc'f/l /ijv(0) = (5 and (R(t)) is the solution to the following SDE 

,22) d «( ( ) = V^W- IT ^ F ffd ( , 

wt/j i?(0) = (B(t)) denoting standard Brownian motion on R and (h(t)) being 
defined in (|20|) . 



Proof. From Relations (|T5|) and (fl7|) . the increasing process of the martingale 

def . / MAr(iV Q *(logiV)t)\ 



is given by 



^fN a * log AT 



(1 + a*)logN Jo a* + 1 + log(Ar- a * + Z N (u))/QogN) 

By PropositionlU this quantity converges in distribution to (2Xt). Hence, using the 
martingale criterion for convergence to Brownian motion (see Ethier and Kurtz |10j ) 
we can conclude that (Mjv(i)) converges in distribution to Brownian motion run 
at speed 2A. 
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With the same notation as in the proof of the previous proposition, the relation 
satisfied by {Z^{t)) yields 

(23) Z N (t) d - (z N (t) - h N {t)) = Z N (0)+M N (t) 

I ]j lcSv { dN ( N ' a * + Zn ^) ~ <MM")) 



(1 + a*) 2 
where 



d N (x) = 



l + (logx)/((a* + l)logN)' 
Let T > and e > 0, and let An be the event 

5~ 



A N = < inf Z N (s) > 

{o<s<T 2 

By Proposition U] and the fact that (h(t)) is nondecreasing with h(0) = S, there 
exists N such that for any N > N , P(A C N ) < e. 
Let us now observe that 

(24 » ^»°H 1+ (a4T)lo sW )}"'- 



On the event An, we obtain from Relation (|23|) and the fact that for N large 
enough, 

; 4 
sup d N < - 

[4/2,oo) ° 

that for any t <T, 

\Z N (t)\ < \Z N (Q)\ + \M N (t)\ + 6{1 l^ )2 J \Z N (u)\du. 
As a consequence, Gronwall's lemma gives us that on the set An, 
sup \Z N (t)\ < ( \Z N (0)\ + sup \M N (t)\ ) e CT , 

0<t<T \ 0<t<T / 

where C — 4/// (2(1 + a*) 2 ). Using the fact that the sequence of processes (Mjv(t)) 
is tight for the topology of uniform convergence, together with our assumption on 
Xjv(0), we obtain that there exists Ni and K > such that for any N > Ni 

P(B c NK )<e, where B NjK = \ inf Z N {s)> 5 -, sup \Z N (t)\ < K 

yo<s<T l 0<t<T 

Next, recall the notation wt for the modulus of continuity of the function /. On 
the event Bn,k, Relation (|23| gives us that for any £ > 0, 

WzJ0<^mJ0 + C sup / [ZNiu^du^wjfJO + CKt 

s,t<T J s 
\a-t\<£ 

The sequence of processes (Mjv(i)) being tight, this relation and the fact that 
P(Bn.k) > 1 — £ show that for any r\ > 0, there exists £o such that for every £ < £o 
and N > N , 
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Since we can apply this reasoning to any e > 0, here again the tightness criterion 
of the modulus of continuity enables us to conclude that the sequence of processes 
(Z N (t)) is tight. 

Let {R{t)) be one of the limiting points. Up to a change of probability space, 
we can assume that the convergence to (R(t)) holds almost surely on [0,T] for the 
uniform norm. But on the set Bn,k, (|24|) and Lebesgue's differentiation theorem 
guarantee that the integral term on the right hand side of Equation (|23p converges 
almost surely to 

Jo K u ) 

Consequently {R{t)) satisfies the SDE (|22[) with R(0) — y, and the uniqueness of 
such a solution gives us the convergence in distribution we were seeking. □ 

A direct consequence of this result is that, starting from JV™ 1 , the process (L^ (t)) 
behaves like an Ornstein-Uhlenbeck process around A" 1 . 

Theorem 1 (A stable Regime for (if (*))). If pi < 1/2, L$(0) = A and Lf (0) 
is such that, for some 

L?(0)-N a " 
lim - = — = y, 

N^ + oo y/N a l log A 

then the sequence of processes 

( Lf (N a "i (log N)t) - N a "i \ 



(25) dZ(t) = ^/2A^dS(^) - j-fl Z(t)dt, Z(0) = y, 



\ y/N a " log A J 

converges in distribution to an Ornstein-Uhlenbeck process (Z(t)), i.e. the solution 
to the SDE 

Mi 

K + i) 2 

where (B(t)) denotes standard Brownian motion. 

To complete this section on the time scale t — > A" 1 log A, the following propo- 
sition shows that, if pi < 1/2 and (0) = A, then the process (L^ (t)) always 
reaches the stable regime around A a i described in the previous theorem. In par- 
ticular, if Lf (0) = 0, the hitting time of [N a * J is smaller than N@ for any j3 > a\. 

Proposition 6. Let 

T N d = inf |s>0: Lf (s) € N a i+ -^N a l logN, N a l logivj X . 
If pi < 1/2, (0) = A and Lf (0) < N? for some /3 G (a>l,l), then 

Proof. As before, the result is proved for the process (X_/v(i)) instead of (L^(t)). 
Suppose that X N (0) = lA* 3 ]. The SDE ([14j) and Relation (fTf]) show that for any 
stopping time t, one has 



(26) E(X N (tAT))=X N {0) 

p / log((l+l K W)r') 



(a* + 1) log A \Jo 1 + a* +log((l + AAr( M ))/A Q *)/logA 



■ dw 
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Defining again 

H x d = inf{s>0:L 1 (s)= L^J 1, 
and setting x q — 2 \N a ] , the above relation gives 

< X N (0) - 1 ° g2 „. -E(tAH x tr). 

~ w (a* + 1) log AT 1 + a* + (log 2)/ log JV v x ° ' 

Consequently, letting t and then N go to infinity yields 

E(J3» 
limsup — 7 , — < +oo. 

N ^+i Nt 3 log TV 

We can therefore assume that X N (0) = 2\N a ""\. Setting this time = [N a ' \ + 
[N a _e J , where e > is such that a* + e < (3, the same argument gives us that 

n < y a io g (i + i/jv e ) , , 

Xl - Xn ^> {a * + i)logNl + a*+ log(l + l/N*)/logN 
and thus 

E(H-n) 

(27) limsup — — 1 < +oo. 

Similarly, if x% = [N a * \ + [N a *~ 2s \ and if we choose X N (Q) = xf , the above 
equation gives for r = H x n 

I jya*-2s I < I AfQ*-e I M log(l + 1/JV 2e ) , , 

L J " L J (a* + l)logAa + a*+log(l + l/A r2E )/logA r V 

so that Relation also holds for H x n. Setting xf = AT* + A^*- ie , we can 
proceed by induction until the smallest integer i* such that i*s > a*/2. Finally, 
we obtain that as A^ —¥ oo, 



E(T N ) HH x n) ^E(H xf -H x 



NP (log A") 2 A^ (log AT) 2 ^ N/ 3 (log AO 



0. 



We conclude by using the Markov inequality. 

Up to now we have been dealing with the case Xjv(0) > A^ Q . There remains to 
consider the case ATjv(O) < A^ Q . First, Proposition [3] shows that we can assume 
directly that Xn(0) = [xN a J for some x € (0,1). We can then proceed as 
before by estimating H Na * _ N c*-e for e sufficiently small and by decreasing the 
exponent by e at each step until it falls below a*/2. This completes the proof of 
Proposition [6l □ 

Remark. 

Proposition [6] completes the results of Section [4] Indeed it shows in particular 
that if if (Q) = and L^{N) = 0, then the average hitting time E(T N ) of the 
neighborhood of A^" 1 is, up to a constant, upper bounded by A"" for any j3 > a\. 
With the same arguments as in the previous proof, it can be shown in fact that 
there exists a constant C > 0, such that 

(28) E (0 , N) (T N ) < Ci bgloglog(JV) ■ 
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6. The Fluid Time Scale 

Recall that the fluid scaling of (L(t)) — (Li(t), L 2 (t)) consists in speeding up 
the time scale of the Markov process proportionally to the norm of its initial state 
and by scaling the state variable by the same quantity. Hence, if \\L (0)|| = 
max(L^(0), L 2 (0)) = N, we are interested in the process 

(L N (t)) d ^±(L(Nt)), 

see for example Bramson j8] and Robert [23] . 

The initial fluid state considered in Propositions |4] and |6] satisfy 

lim L N (0) = (0,1). 

N— >+oo 

It has been shown in Proposition [5] of Section [5] that, the hitting time of N" 1 by 
Li(t) is only 'slightly' larger than 7V Ql . The following proposition completes this 
result on the fluid time scale. 

Proposition 7. Suppose that p\ < 1/2, p 2 > 1/2, and that (L^(t),L2 (t)) is the 
solution to the SDE (|12[) with initial conditions (0) = N and (0) such that 



lim L?(0)/N = x€ 

N— >-+oo 



Let Tn be defined by 



T N =miis>0: Lf (s)€L 2 v (s) a 



-y/L»(s)°l logics), yfl$(Jr l logics) 



Then, for the convergence in distribution we have 

T N def. 2x 



n^+oo N pi — 2Ai 

Proof. To start with, note that for both i £ {1,2} and all t > 0, L0) < L,(0) + 
■^VxiCfO, Hence, by the law of large numbers for Poisson processes, for every 
j] > and K > there exists Nq £ N, such that with probability greater than 
1 — 77, the relations 

Li(tN) <(x + 2\ 1 K)N and L 2 (tN) < (1 + 2\ 2 K)N 

hold for all t < K and N > N . 

Let us now define = inf{s > : Lf (s) < AT/ (log iV) 2 }. Of course, this time 
is when L± (0) < N/(\ogN) 2 . By the remark made in the previous paragraph, 
with probability at least 1 — r\ we have for every s < A (KN) 

logjV-21oglogiV logLfCg) 

> ' 1 /1io\ J&"\ oi„„i — at i 01 — i\r — 1 r JV I „\ 1 1„„ r iV 



and 



log(l + 2A 2 A:)-21oglogiV + 21og7V " log if (s) + logics) 
logLf(s) log(l + 2A 2 A-)+log^ 



logif(s) + \ogLg( s ) " log(l + 2A 2 A-)-21oglogA^ + 21og7V' 
Hence, on this time interval the service capacity offered to class 1 jobs (resp., class 2 
jobs) is at least (resp., at most) 1/2 in the limit. Since p 2 > 1/2, this implies that 
the process (L 2 (s)) is increasing on the time interval [0,^ A (KN)]. That is, for 
N sufficiently large we have with probability greater than 1 — r\ 

inf LF(s) >N- log N. 

s<Tf,A(KN) 
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Consequently, Relations (|29|) can be completed by the inequality 
log if (s) log(i + 2X 1 K) + log AT 



< 



logif(s) + log if (s) " log(a: + 2AiJQ + log(l - (log TV) /TV) + 2 log A^' 

This shows that, with high probability, as A~ gets large the two queues receive 
the capacity 1/2. As a straightforward consequence, we have the convergence in 
distribution 

lim %- = ^— = t (x). 

jv^oo AT nx - 2Ai v ; 

Next, let us suppose that if (0) < |A7(logA0 2 J and let us define 

= inf {s > : if (s) < 2N a i}, 

Since if (A*i) grows at linear rate (recall that pi > 1/2), we can again compare 
if to A'jy and conclude from Relation (f2l)l) applied to the stopping time Tjy that 

Ml l0g(2) E(r^)<if(0). 



(a* + 1) log N 1 + a* + log(2)/ log AT 
Consequently, 

lim E (rir) /A" = 0. 

JV->+oo v 7 

We can thus assume that if (0) = 2[A ra iJ, and Proposition [6] shows that in this 
case, 

lim E{T N )/N = 0. 

N— >oo 

Coming back to the initial question (with an arbitrary x) and using the strong 
Markov property of L N combined with the last two limits, we obtain that T/v = 
r N + (rf - r«f ) + (T N - rf ), where 

— ^*o(z), — >0 and — >0 

as N co. By Slutzky's Lemma, we can conclude that T/v — >■ to(x) in law. □ 

Theorem 2. Suppose that p\ < 1/2, pi + p2 < 1 arcd iei k > 0. Then, there exists 
t]q > suc/i i/iat /or any sequence (Zf ) satisfying 

jN 



lim sup 

W-S- + 0O 



- 1 



ak 

i/ (if (t), if (t)) is the solution to the stochastic differential equation (|12[) with 
initial condition (ii(0), i2(0)) = (l^,N), then 

if( S ) 



(30) lim P sup 



- 1 



> K = 0. 



(if( S ))^ 

Proof. Let us write 5 = 1 + k, and let us prove that there exists 770 > such that 

N~+oc" \ <s<7 10 N (if (s)) Q J 



(31) Jim P ( sup /TN \\i, >S)=0. 



The other inequality can then be shown by using the same technique, and so we 
omit the details. 

/y 

The process (if (s)) is stochastically bounded from above by (i 2 (s)) describing 
the number of class 2 jobs in the priority system where class 1 jobs have the priority 
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of service, i.e. queue 2 is served only when queue 1 is empty. For this system, it is 
not difficult to show that the convergence in distribution 

lim P 2 8<l)=(l + M2(pi +P2- l)s, s < 1) 

holds. Similarly, the process (if (s)) is stochastically bounded from below by 
(if (s)) describing the number of class 2 jobs when they have priority, and one has 
the corresponding convergence in distribution 

(L»{Ns) \ 
J iT t = ^7 — , s<l \ =(l + fx 2 {p 2 -l)s, s<l). 

Fix < rjo < I such that 

5(1 + 2 M2 (p 2 -f)7 7o ) Q i > f. 
Let also £ be the event defined by 

. f LP (a) „ if(s) 3 , 

Lo<s<i7 A r ^ V 0<s<jj o JV iV 2 

The above convergence in distribution results show in particular that, for any e > 0, 
there exists No such that for every N > Nq, 

V(£C)<p( sup hM >1 \ +P ( M ffi!l<i+3 (pa _i)J < e . 
^o<s<i)ojv ; \o<s<J7oiv 7v 2 y 

By Relation p7|) . the service rate of queue 1 at time s is given by 
def. x log[(l + L?(s))/(L? ( S ))<] 



A(a) = Ai+/ii 



+ f )((«; + 1) log(if (*)) + log[(l + if ( S ))/(if (s))^]) ' 

If for some y > 1 and some s < 1 if (s) > y(if (s))" 1 , then 

An > A log(y) 

U_ 1 Ml K + l)((aJ + l)log(Lr(s))+log(y))" 

Furthermore, 

(32) A( S )>M2/) d = Ax+W- l0g(y) 



K + l)(K + l)log(iV)+log(y)) 

holds on the event £. 

Denote by (X y (u)) the birth and death process on Z whose +1 [resp. —1] jumps 
have rate Ai [resp. pn(v)} and starting at 0. The corresponding reflected process 
is denoted by (X+(u)). This is in fact an M/M/l queue with input rate Ai and 
service rate /Ltjv(y)- 

Define 

oi = o\l + -P2{P2 - Ij^oJ and 8 2 = — - — . 

Since p 2 < 1, the definition of 770 gives that 81 > I and therefore (52 > 1- Note that 
8 > 8 2 . Suppose first that Zf /iV Q i -> (5. 

A simple coupling argument gives that the processes (if (s)) and (X<s 2 (s)) can 
be constructed so that on the event £ , the relation 

(33) if + 

holds for every s < inf{u : if (u) / (L% (u)) a * < 8 2 }. 
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For every x > 0, Kingman's inequality (see again Kingman |15j or Relation (3.3) 
of Theorem 3.5 in Robert [23]) gives us the estimate 

(34) P ^supX & («) > x^j < exp(-x(n N (5 2 ) - A x )) . 

In particular, the random variable 

^sup^W 

converges in distribution to since 

Mi log(<5 2 ) 



UN (62) ~ Ai + 



{a\ + l) 2 logA^' 

as N goes to infinity. Let 

T N = inf {s > : if + X 5a (s) < (5 2 iV Q i}. 

In a similar way as in the proof of Proposition [B] for example, a simple drift analysis 
shows that 

vlT s . K + i)(K + i)iog(AQ + iog(a2)) , N 

Mi log(oa) 



and consequently 



lim ¥(T N > AT/ (log AT)) = 0. 

./V— >+oo 



From Relation (|3"3"|) . we can conclude that with probability tending to 1, if /N ai 
does not grow above if /N ai until the time T/v, which itself occurs much before 
N. As a consequence, it is enough to prove Identity (|3~Tj) with the assumption that 
Lf (0) = AT and 

hm ^°>=fc. 

As before, thanks to Inequality (1521 . one can construct a coupling such that the 
relation 

(35) Lf(s)<Lf(0)+X+( S ), 

holds for every s < 770 A^ on the event £. 
For every y > 0, let us define 

f L N (s) 1 

r, = inf{ S >0: X+(*)>y} and fT = inf | s > : ^LAJ_ > ,5 J. 

The proposition will be proved if one shows that ¥(H < r/oN) converges to as N 
goes to infinity. 

On the event £, if < s < tjqN is such that Lf (s) > (5(if (s)) a i then 

> <5(l + \^{pi ~ l)vo) ai N a i = 5iN a '. 
Consequently, Equation (|35|) and the definition of 82 give that for every N > No, 

(36) ¥(H < Vo N) <e + P(r r ; (ji _ 1)/2 < %iV). 

But using the same technique as in the first part of the proof of Proposition [4j we 
can show that the probability that reaches iV" 1 (5\ — 1)/2 in one excursion away 
from is so low, that the probability that it reaches this height during one of the 
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O(N) excursions is does in the interval [0,rjoN] tends to as N tends to infinity. 
Therefore, we can conclude that 

F(H < T} N) < 2s, 

for TV large enough. Since this property holds for every s > 0, Theorem [2] is 
proved. □ 

Corollary 1. Under the assumptions of Theorem^ the convergence in distribution 



lim 



N 



< t < t = ( 7 (t), < t < to) 



holds, with t — l/(/i2(l— p\ — pi)) and ^(t) = 1 + P2(pi + P2 — !)£■ 
In addition, for every t < to we have 

L?(Ns) 



(37) 



lim 1 

7V->+oo 



sup 

0<s<t 



- 1 



> K = 0. 



Proof. Let us first prove the convergence 

(38) lim ( L "[ Nt) ,0<t< m j = r.i/i.l)-. / ; ,/,, i 



N 



The SDE (| 12[) is used in the same way as before, and so we only sketch the proof. 
By Theorem [21 we have 

L?(Na) 



lim 1 



sup 

0<s<r;o 



- 1 



< K 



1. 



Thus, Lf(Ns) is at most of the order of N" 1 with arbitrarily large probability. 
This implies that for any s < 770, all the arrivals at queue 1 up to time Ns are 
processed. Hence, queue 1 uses the fraction pi of the capacity of the server, and 
the remaining capacity is devoted to queue 2. The convergence Q38p is proved. 
Furthermore, from Relation (1301) we obtain that 

L?(N Vo ) 



lim 1 

N-t+oo 



(L«(N Vo )ri 



- 1 



< K = 1, 



and L2 (Nijo) ~ j(r)o)N. Consequently, Theorem [2] applied with the initial condi- 
tion (L?(Nrio),Lg(Nr)o)) shows that the convergence Q38|) and Relation p0|) can 
be extended to the interval [0, 770(1 +l{vo))]- That is, 



(39) 



lim I 

7V->+oo 



sup 

,0< S <7jo(l+7(>)o)) 



Lf (Ns) 



(L?(Ns))<*l 



> K = 0. 



Proceeding by induction, as long as 7(2;) 7^ if the convergence (|38| and Inequal- 
ity ([39]) hold on [0,x], these relations can be extended to [0,x + 7(^)770] . The 
corollary is proved. □ 

The following theorem is the main result of this section. It shows that, on the 
fluid time scale, the process (L^(t)) lives on the space scale N" 1 as long as the 
process (t)) has not reached 0. 

Theorem 3. Suppose that p\ < 1/2 and pi + P2 < 1. If (-^1 (*)>-^2 W) * s the 
solution to the SDE (12} with L%(0) = N and 



Lf(0) 



lim 
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then we have for the convergence in distribution 
~L?(Nt) Lj?(Nt) 



lim 



0<t<t o ) = ([ 7 (i) ai ,7W],0<t<io), 



N a l ' N 

where to = 1 / '(^(l- P\ — Pi)) and j(t) = 1 + p2(pi + P2 - l)*- 

Proof. From Corollary[TJ we obtain that Relation (1571) holds any k > 0. Hence, the 
process 

WWW*' °- s<t ° 

converges in distribution to the process constant equal to 1 on [0,to) as N goes to 
infinity. We conclude thanks to the convergence of (L% (Nt)/N, < t < to). □ 

We can now state a fluid limit result concerning this network under the assump- 
tion pi < 1/2 and p 2 > 1/2 and with the more general initial conditions considered 
in Section SJ Analogous statements for the other cases are available, but for the 
sake of simplicity we do not give them here. 

Corollary 2 (Fluid Limits). Suppose that pi < 1/2 and p 2 > 1/2. If(L?(t),L$(t)) 
is the solution to the SDE (11211 with initial conditions such that 



lim 

7V-> + oo 



Lf(0) 



x and lim 

N Af-^+oo N 

for some x € [0,1], then for the convergence in distribution we have 



^(0) 



1 — x, 



lim 



(*)), 



TV ' N 

where the pair (^1,^2) is defined as follows: if to(x) = 2x/(fii—2Xi) 



Ar 



Pi 



t, (l-x)+ 



A 2 



P2 

2 



0,(l 2 (t )+[A 2 -M2(l-pi)] (t-<o)) H 



t>t . 



Proof. We give only a sketch of the proof. Until the time Nto, both queues are of the 
same asymptotic order for the log function. Consequently, they both receive half of 
the capacity. Notice that on the time interval [0, Nto], because of our assumptions 
the variable L± decreases whereas L 2 increases. After time Nto, the variable 
is of the order of N a i and is therefore negligible for the fluid scaling. □ 

Remark. 

The constant a* does not show up explicitly in the fluid limit, but clearly enough 
this result is much weaker that Theorem|3] As we shall see in Section[7J the constant 
a* plays an important qualitative role in the expression of the invariant distribution 
when the regime is close to saturation. Hence, a\ is not an artificial artifact of the 
fancy time scales used in this analysis but is definitely an important feature of the 
network. 
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7. Heavy Traffic Regime 

When pi + p-2 < 1, since the queueing system is work conserving the Markov 
process (Li(t), i 2 (t)) has an invariant distribution 7r p , where p = (px,p%). I n the 
following, (L Pi i,Lp 2) stands for a random variable with distribution ir p . Obtain- 
ing explicit expressions to describe n p seems to be quite difficult: because of the 
logarithmic weights, the double generating function of (n p (to, n) , (to, n) € N) does 
not satisfy an autonomous equation as it is often the case, e.g., in classic two- 
dimensional reflected random walks. A precise result on the asymptotic behavior 
of tt p when p\ + P2 is close to 1 can nevertheless be obtained. Its proof relies heavily 
on the scaling results of the previous sections. 

Theorem 4 (Heavy Traffic Limit of Invariant Distribution). If p = (p\,p2) and 
pi < 1/2 is fixed, then as p — pi + P2 tends to 1 we have the convergence in 
distribution 



lim 



1 ((l-p)^L pA ,(l-p)L p , 2 ) ^(Ef\E v ), 
1-P1) ^ ' 



where here again a* = p\/{l — pi) and E v denotes an exponentially distributed 
random variable whose parameter n is given by 




Proof. For every £ > 0, let E^ denote, as before, an exponential random variable 
with parameter £. All the variables used here will be assumed to be independent. 
The workload of the system has the same (invariant) distribution as the workload 
of an M/G/l queue, with arrival rate Ai + A 2 and with the same service distribution 
as 

<r*£ BE^+(l-B)E P2 , 

where B is a Bernoulli random variable with parameter Ai/(Ai + A2). Hence, 
Kingman's heavy traffic result for the workload at equilibrium (see Kingman [16] 
or Proposition 3.10 of Robert [53]) gives the convergence in distribution 

(40) lira (i-p)(bhL + hlA =E ^ 

A 2 /V 2 (l-Pi) V Ml M2 / 

where 770 is the constant given by 

Vo= lim ; 2 =V2 / Jl+ { ^-^ )2 Pl {l- Pl ). 

A2/V2U-P1) (Ai + A 2 )vw - #a 1+ a 2 ) / y P1P2 

In particular, the family of random variables 

(41) [(l-p)(L Pil ),(l-p)(L p , 2 )] 

is tight as p /*• 1. Let us denote a possible limit by (Xi, X2), corresponding to a 
sequence (p n ) converging to 1. For every n > 1, let us define (L Pnt i(t), L Pn ^(t)) as 
the process with initial condition (L Pnt x,L Pnt 2), which is thus a stationary process. 

The strategy of the proof can be described as follows: use the results of Sections[4] 
and[6]to prove that, for a fixed time t n , L Pn ^(t n ) is large and £ P „,i(£n) is close to 
(L Prei2 (^ri)) Ql ■ Consequently, by stationarity the equivalence L Pnt \ ~ L a p ^ 2 holds 
and Relation (|4H|) then gives us the desired result. 
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To start with, let us assume that P(X 2 = 0) > 0. Then, there exists 6 > with 
the property that for any rj > 0, there exists no £ N such that for every n > no, 

(42) P fi P „,2 < tz^-) ^ S - 

V 1 Pn J 

Now, Relation (gO) tells us that P(Xi + X 2 > 0) = 1. Consequently, for any e > 
there exists 770 > and no such that if n > no and 771 + 7/2 < 2770, 

(43) P < 7^,^,2 < j^-) < e. 

V 1 - 1 - Pn J 

Let so > and 771 > 0, and set t n = so/(l — /O n ). For 71 large enough we have 

(44) 5 < P fi p „ )2 < T ^-) = P (L Pnt2 (t n ) < 

< P ( L Pn , 2 (t n ) < — sup L Pn , 2 (s) > , n 770 _ J +£ 

\ 1 - Pn [o,t„] ^(l - p„j y 

+ P (L p „, 2 (i n ) < sup L Pni2 (s) < Vo _ - y L Pn<1 (0) > j . 

\ 1 - Pn [0,i„] 3(1 - p„) t~PnJ 

Clearly enough since (L 2 (i)) decreases at rate at most /j 2 , as n increases the first 
term on the right hand side of (|44| can be made arbitrarily small by choosing 771 
and so m such a way that 

(Ci) 771 < y + (A 2 - /j 2 )s = y - /J 2 (l - p 2 )s . 

The third term in the right-hand side of Relation (|44[) can be upper bounded by 

(45) K L ^ 0)> T^'& L ^ a)< W^. 

(46) Zw(0) > inf L Pnil (s) > 2% _ 

I- p n [0,t„] 3(1 - p„) 

As before, the quantity in (|4"5|) will tend to as n — > 00 if we choose so such that 
(Ca) j-Mi(1-Pi)*o>0. 
Finally, if 

sup i p „, 2 (s) < 770 _ and inf i p „,i(s) > r ; 
[o,t„] 3(1 - p n ) [o,t n ] 3(1 -p„) 

the infinitesimal drift A Pni2 (s) of £p n>2 satisfies on [0, £„] 

A M> f fog (770/(3(1 ~ P»))) \ / 1\ 

Ap„,2(.) > P2 ^P2 - log(W( 3(l-^))) + log(2V(3(l-p„ )) )J ~ ^ v ^ 2 - 2 J ' 

as 77 goes to infinity. But p n>2 converges to 1 — p\ > 1/2 by assumption, which 
means that the quantity in (l4l)|) will tend to as n gets large whenever 

(C 3 ) 771 < /i 2 (l/2-pi)s . 
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Choosing first so small enough to match (Ci) and (C2), and then rji small enough 
to satisfy C\ and C3 , we can conclude that the right hand side of Relation (|44"1) can 
be made smaller than <5/2, which yields a contradiction. We have thus proved that 
¥{X 2 > 0) = 1. 

As a consequence, for any e > 0, one can find 770 and no G N such that for every 
n > no, 

The last step follows an analogous path, in that we choose a convenient time t n to 
first show that, for any /3 > a*, the quantity 

T > 7]1 T > 1 

1 - pn (1 - Pn) P 

is arbitrarily small with the help of Proposition [BJ and next that for K sufficiently 
large, 



We (Lp„, 2 ) Ql 



-K a/ (L p „, 2 p log (L Pn , 2 ), KJ (L p „, 2 ) Ul log (L Pn , 2 ) 



with high probability by Proposition [7] 

Since L Pnt i is of the order of {L Pn ^) ai and that a\ < 1, Relation pt)|) gives the 
desired convergence result. □ 

8. General Functions for Resource Sharing 

In this section, we consider the case where the function (log a;) describing the 
access to the resource is replaced by some other function (f(x)). For a two node 
network, the corresponding Q-matrix is given by: for every x e 

!q(x, x + = Xi 

q(x,x-ei) = fit 



f(xi) + f(x 2 y 

To concentrate on the most interesting case, throughout this section we assume 
that pi < 1/2. We analyze only the first two time scales in order to stress the 
difference with the log function, at least concerning the scaling parameters. The 
proofs of the results presented here are similar to those in the log case, and so most 
of the time we only sketch them. 

8.1. The log log Function. Here we consider the function f(x) = log(log(e + x)). 

To simplify the notations, we use instead the function x —t log 2 (x) d =' log(loga;). 
This, of course, does not change the convergence results obtained in the following 
paragraphs. 

Initial Phase. The first time scale is tt-$- cf>jsf(t) with 

^(t) = exp[(logiV) t ] . 
The stochastic evolution equation is in this case 

(48) Li(<M*)) 

= Lj(0) + AxM*) " Mi [ Mt) ^ 2 ^ff } du + M N (4> N (t)), 
Jo log 2 (Li(u)) + log 2 (iV) 
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where (Mjv(i)) is a local martingale. Since (f>N(t) <C iV as long as i < 1, the second 
coordinate L2 stays at N at least up to t s» 1. This justifies the fact that in the 
above expression, the term log 2 (L,2(u)) has been replaced by \og 2 (N). 
Let us now define Z^{t) = £1 (0at (£))/<?!> at (f). We have 



Z A r(t) = Z A r(l) + A 1 - 

Mi 



1 



<M*) 
t 



\0g 2 (Z N (u)(j)N(u)) 



(j> N {t) Jo \og 2 {Z N {u)4> N {u)) + \og 2 (N) 
Note that, for every u > 0, 



$v(u) dw + 



M N (^ N {t)) 

<M*) 



\og 2 (Z N (u)cj> N (u)) = \og 2 (<j> N (u)) + log 1 



log(Z jv ( U )) 
log(0jv(u)) 



and 



log 2 (0jv(w)) = ulog 2 (iV). 

Hence, using the same methods as in Sectional we obtain the following equivalence 
for the convergence in distribution of processes, uniformly over any compact subset 
<rf(Q,a;): 



(Z N (t)) ~ Ai 



Mi 



(*) 



— 0^r(w) dw 



A, 



-Mi 



t+1 



The following proposition is the analogue of Proposition [5] for the log log function. 
Recall that a* = pi/(l — pi) < 1 since p\ < 1/2. 

Proposition 8. If (L?(t),L$(t)) is the Markov process with Q -matrix (|4T[) and 
with initial condition (0,7V), then the convergence in distribution 



lim 



Lf (exp [(log Ap *]) 
exp [(log A)*] 



, = (Ai - Mi — -j-jO < t < a\ 



holds for the uniform topology on compact sets of (0,0^). 

The Local Equilibrium. The last proposition together with the results obtained in 
Section [5] suggest that, if p\ < 1/2, the process should remain stable around the 
value exp [(logiV)" 1 ] . As in Section El let us assume that Lf(0) — S<j)N(a*) for 
some 5 < 1, while L 2 (0) = N. Using the relation p\ — a*/(l + a*) and the 



■ du + M N (t). 



evolution equation B8|) . we obtain that for every t > 



Lf(t) = Lf(0) 



A' / 



Mi 



log[log(if ( U ))/(log7V) Q i 



1 y log[log(Lf (u))/(logJV)«I] + (at + l)log 2 (A) 

As we shall see, the appropriate scaling of time around ^(a*) turns out to be 
1 1 y ipNt, where 



Vjv = </»ivK)(logiV) a Mog 2 (iV) 



exp 



(logA) Q i (log7V) Q Moglog(7V). 



Indeed, let Zjv(t) = (ipNt) / '4>n (a*) ■ For every it > 0, we have 



log 



log(Lf(V>Am)) 
(logA) Q I 



log 



\og(Z N (u)) 
(log A) Q ? 
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so that 

M N (i; N t) M N (1) 



Z N (t)-Z N (0) . 

4>N{al) <PN(at) 

Mi f^ Nt log[log(Lf( U ))/(log7V)^] 



(at + 1)<M<**) Jo log[log(if (u))/(logN)<*l] + (a\ + l)log 2 (iV) 
Ml (log7V) Q Mog 2 (iV) /■* log^vH) 



K + l) Jo log(Z N (u)) + (al + l)(logN) ai \og 2 (N) 



du. 



With the same tightness argument as in the proof of Proposition HJ we obtain the 
corresponding convergence result detailed below. It is remarkable that the limit is 
the same as in Proposition |4j 

Proposition 9. If (Lf (t) , L$ (*)) is the Markov process with Q-matrix (|47p . and 
initial conditions L 2 (0) = N and L± (0) ~ <5exp[(logiV) Q i] for some S £ (0,1], 
then for the convergence in distribution of processes we have 



N- 

where 



lim 



n (if (^ N t)e-^ at ) = (h(t)), 



(logJV)^l (logiY)" 1 log log N 



-0AT = exp 



and (h(t)) is the function defined by Relation ([20 

-N 



Furthermore, if L± (0) ~ exp[(logiV) Ql ] + y^Ji^N for some y € R, then th 
sequence of processes 

( L»ty N t)-eV°* N ) at \ 
{ ) 

converges in distribution to the Ornstein-Uhlenbeck process defined by Relation (|25 



e 



Remark. 

If p\ < 1/2 and if the initial condition is (0, JV), we obtain that L\ is of the order 
of 

exp[(logiV) Q i] 

on the time scale t i-> tpNt- This is much smaller that the quantity JV" 1 corre- 
sponding to the log policy. One can then take a function / growing more slowly to 
infinity, such as log log log. Under the same assumptions, the variable live in 
a region with an even smaller order of magnitude. By pushing this scheme a little 
further, we would obtain a policy similar to the Head of the Line Processor-Sharing 
(see Bramson [7]), where node j receives the bandwidth 



l{m#o} H 1" 1 {n J ^o] 



The main drawback of this policy is that a node with many jobs has the same 
fraction of the capacity as a node with only a few of them. For this reason, it is 
more likely that local congestion will occur more frequently 
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8.2. The Time Scales for a General Function. Finally, let us return to the 
general case. Let us assume that the function x i— > f(x) is an increasing continuous 
function on R + , tending to infinity as x — ¥ oo, and suppose that there exist two 
functions x i— > Af(x) and x t— > Bf{x) on R + such that for any t > and z > 0, 

(Fl) hm = 0, hm — ^ = Btiz). 

The first time scale for the initial phase is given by 
(49) <j> N :t^ f-\tf(N)), 

where f" 1 denotes the inverse function of /. If (L^ (t), L% (t)) is the Markov process 
with Q-matrix (|47|1 and initial condition (0, N), then the convergence in distribution 

Hm / L? {ct>M{t )) >Q<t<a *) = ( Xl _ ^ * Q < t < a * 



N^ + oo \ (j) N {t) ' V V * + l 

holds for the uniform topology on compact sets of (0, a*). Indeed, the first relation 
in Condition (F\) ensures that the second coordinate stays at N while is 
of the order of 0Ar(i), so that exactly the same techniques as in Section [4] can be 
applied. 

The second time scale of interest is 1 1— > ipNt, where 

50 = . 7-7 7 ^ ■ 

If [Lf{t),L^{t)) is the Markov process with Q-matrix (|47|) and initial conditions 
£2(0) = N and L\{Q) ~ S4>n(<^i) for some S £ (0, 1], then we conjecture that the 
convergence in distribution of processes 

n^+oo \ cl> N (a$) J V V >h 
should hold, where (h(t)) is the function defined by 



(51) 



h = 1 if 6 = 1, 

h(t) 1 ut 

7~T dlt = — -rr if 5 7^ 1 . 

B/Cu) (1 + ^) 2 



Some regularity properties (required to use Lebesgue's Differentiation Theorem for 
instance) are clearly necessary to justify this convergence, but then the rest of the 
proof should follow the lines of the proof of the corresponding result for the log 
case. 



Some Examples. 

a) f(x) = loglog(ir). 

In this case Af(x) = 1/logx and Bf(x) — logs, and we recover the expres- 
sion of the time scales obtain in Section 18.11 for the log log function, that 
is 

cf, N (t) = exp[(logA0*] and i/j n = (\ogN) a i loglog(JV) exp ((log7V) Q i) . 
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b) f(x) = (log a:)". 

Here Af(x) — l/[/3(logx)' 3_1 ] and Bf(x) = \ogx, we gives the two time 
scales 

4> N (t) = N' 1 " 3 and Vjv = paf-^N^QogN)^- 1 . 

In particular, we recover the time scales of the log case which have already 
been identified. 

Note that the functions x h-> x a , a > 0, do not satisfy the first relation in Condi- 
tion (Fx). 

9. The Network with J Nodes 

In this section, we briefly describe the case of J + 1 nodes competing for the 
single resource. A special case is considered to illustrate the similarities and also the 
differences in qualitative behaviors. The motivation of this section is to show that 
the analysis of the two node network gives the main ideas to start the investigation 
of more complicated situations. New difficulties are indicated in the text, and the 
proofs of these results will be the subject of a further work. 

Let us assume that 

.7+1 

pi < p 2 < ■ ■ ■ < pj+i and ^ Pi < l - 

l 

If the state of the system is {nj)i<j<j+i, the fcth station receives the fraction of 
service 

log(l + n k ) 

log(l + m) + log(l + n 2 ) + ■ • • + log(l + nj+i) ' 

From now on, we consider the case where Lj (0) = for j = 1, . . . , J, and (0) = 
AT. As before (L^t j) denotes the Markov process describing the number of requests 
in each queue and with this initial condition. 

9.1. Initial Phase. The following Proposition is analogous to Proposition^ and 
can be proved in the same way. 

Proposition 10. If [L^ (t)) is the solution of the SDE (fT2)) with the initial condi- 
tion Lj(0) = for 1 < j < J, and Lj+i = N, and if 

def. Pi . 1 

h = I" Jp- X < 1; l - e - Pl < 7TT' 

then, for every 1 < j < J, the convergence 

'L?(N*) \ ( t 
lim I - J —— — , < t < ti = A 7 - - Uj- , < t < ti 

N^+oo \ N* J V Jt+1 

holds for the uniform topology on compact sets of (0,ii) 



Observe that the quantity t\ is precisely the first time t for which 
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9.2. Second Phase. Let us now give a more heuristic description of the evolution 
of the network after "time" TV* 1 , but still on the time scale i i— > TV*. As we shall see, 
for the second phase the exponent in N of the random variables (Lj(N t ), 2 < j < J) 
are still t, but the exponent 0:2,1 (i) °f ^1 becomes a linear function of t with slope 
less than 1. 

Let us use the different phenomena observed in the two node case to infer the 
behavior of the (J + l)-node system. 

First, since pj > pi for every j > 2, the infinitesimal drift of Lj(N t ) remains 
positive at least for a small amount of time after t\. Hence, one should have the 
following convergence in distribution: for 2 < j < J, 



lim '■ 



Lf(JV*) 



TV* 



,h<t<t 2 



A, 



Pj 



t 



a lt2 (t) + (J - l)t + 1 



ti < t < t 2 



where t 2 is the first time t at which 

A 2 = P2 



t 



a lt2 (t) + (J - l)t + 1' 

Second, the station 1 should remain at a local equilibrium in the sense that the 
coefficient ai j2 (t) should be determined as follows: 

1 ^ a h2 (t) + (J - l)t + 1 
Combining these relations, we obtain that 



ai, 2 (t) = 



Pi 



-(1 + (J - l)t) and t 2 



P2 



1-Pl 

Of course, t 2 has to be strictly less than 1. 



l-pi-(J-l)p 2 ' 




Figure 2. The network with J + 1 nodes on the t ^ N l time 
scale when tj < 1. 



9.3. Subsequent Phases. Let us now give a, still heuristic, description of the fcth 
phase for 1 < k < J. The first k — 1 stations are at a "local" equilibrium. Denoting 
the exponent of the jth station in the fcth phase by Oj^t), the equilibrium is 
characterized by the relation 

Qj.fcft) = 

<*!,*(*) + a 2 , k (t) + ■■■ + a k -i, k (t) + (J — k + l)t+l Pj ' 
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lox 1 < j <k — l. The end of the fcth phase, tk, corresponds to the situation when 
the fcth station reaches an equilibrium, i.e. tk satisfies 

tk 

o:i,fe(*fc) + ot2,k(tk) H 1- ock-x,k{tk) + {J - k + l)t k + 1 

Consequently, we obtain that 

«i,fc(t) = 5^(1 + (J-k + l)t), 1 < j < k - 1, 

1 - £i=l P^ 

and 

. Pk 
tk 



1 - Eti Pi-(J~k + l)p k ' 



The time tk is such that tk < 1 if 

fc-i 

^Pt + (J-k + 2)p k < 1. 

9.4. Final Phase. Provided that tj is strictly less that 1, at the time N tj all the 
stations are at a local equilibrium around N aj ' J+1 where i is given by 

Pi 

1 - E i= i Pi 

Keep in mind that strictly speaking, the local equilibrium of the jth station is on 
the time scale t H> N aj - J+1 log Nt. Hence, the whole process can be thought as a 
collection of stationary processes evolving on different time scales. See Figured] 

The difficult technical problem to solve here is for the J—l phases in between, 
during which some of the exponents depend on time, adapting to the linear growth 
of the other exponents. One of our main problems throughout this work has been 
that we did not succeed in proving convergence of the log(Lj)/ \ogN variables 
without showing a more demanding result, namely a convergence result for the 
variables Lj. This difficulty is even more serious here since the exponents depend 
on time. 
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